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We study the properties of compacton-anticompacton collision processes. We compare and con- 
trast results for the case of compacton-anticompacton solutions of the K{l,p) Rosenau-Hyman (RH) 
equation for I = p = 2, with compacton-anticompacton solutions of the L{l,p) Cooper-Shepard- 
Sodano (CSS) equation for p — 1 and I = 3. This study is performed using a Fade discretization 
of the RH and CSS equations. We find a significant difference in the behavior of compacton- 
anticompacton scattering. For the CSS equation, the scattering can be interpreted as "annihila- 
tion" as the wake left behind dissolves over time. In the RH equation, the numerical evidence is 
that multiple shocks form after the collision which eventually lead to "blowup" of the resulting 
waveform. 
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I. INTRODUCTION 

Manifestations of traveling-wave solutions with com- 
pact support arc abundant in nature [i- ]. Nonlinear 
phenomena with linear dispersion lead to solitary waves 
(solitons) that are highly localized in space and have 
wings or tails that extend to infinity. In turn, the interac- 
tion between nonlinearity and nonlinear dispersion gen- 
erates exact compact structures (compactons) that are 
solitary waves defined on a compact support, i.e. soli- 
tons free of exponential tails. Both solitons and com- 
pactons have the remarkable property that they scatter 
elastically, that is, they reemerge with the same coherent 
shape after a collision. However, the scattering of com- 
pactons also gives rise to a compact ripple. The positive- 
and negative-amplitude parts of the ripple decay slowly 
into low-amplitude compactons and anticompactons, re- 
spectively. 

Rosenau and Hyman (RH) discoverred the compactons 
in their studies of pattern formation in liquid drops us- 
ing a family of fully nonlinear Korteweg-de Vries (KdV) 
equations [ ], 



0, 



(1) 



where u = u{x,t) is the wave amplitude, x is the spa- 
tial coordinate and t is time. Following their original 
discovery, compactons have been the subject of intense 
scrutiny [8-24]. A recent review of nonlinear evolution 
equations with cosine/sine compacton solutions can be 
found in Ref. 

Because the K{l,p) compacton equation, Eq (1), does 
not exhibit the usual energy conservation law, F. Cooper 
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and his collaborators proposed a series of Lagrangian gen- 
eralizations of the KdV equation [26-32]. The first such 
generalization was proposed by Cooper, Shepard and So- 
dano (CSS) based on the first-order Lagrangian [2G] 



£(/,p) = 



dx, (2) 



which leads to the CSS equation, 
1 



L{l,p) : ut + 
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with u = (px- We note that in the RH equation, Eq (1), 
K{l,p) corresponds to the set L{1 — l,p+ 1) for the CSS 
equation and in particular K{2, 2) solutions in the RH 
equation correspond to £(3, 1) solutions in the CSS equa- 
tion. For I — p the RH compactons have the property 
that the width is independent of the amplitude. Note 
that {I = 2,p = 1) is the KdV equation. 

In general, for the RH equations there are two conser- 
vation laws. 



Ju{x,t)dx, Q = J uP+^{x,t)dx . 
The CSS equation has three conserved quantities 
AI — I u(x,t)dx, P— I u^{x,t)dx, 



(4) 



E = 



u{x, t) dx , 



1(1-1) 



- auP{uxy 



dx . 



(5) 



To study the solitary waves, one lets u — f{y), where 
y = X ~ ct, in Eq. (1) to obtain (/' = df /dy) 



dy 



d^ 
dy 



(6) 
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Integrating once we obtain 



(7) 



Multiplying by dfP/dy and integrating we obtain the 
equation 



P+l p+l 2 



^{F)Y+CiF + C'2, (8) 
dy J 



where Ci and C2 are integration constants. The com- 
pactons are solutions with Ci = C2 — 0; f satisfies the 
differential equation 



2c 



f3-p_ 



(9) 



p{p + 1) p{p + /) 

For / = p < 3 the solutions are 

f ^ Acos^'^^-^\l3{x~ct)], -7r/2 < /3y < 7r/2 , (10) 
and zero elsewhere. In particular for — (2,2) 



/3 = l/4, A=-c. 



(11) 



The compactons are weak solutions being a combi- 
nation of a compact object f{x) confined to a region 
(say initially —xq < x < xq and zero elsewhere). At 
the boundaries ixo the function is assumed continuous, 
/(xq) = 0, but the derivatives most likely are not. For 
there to be a weak solution we need that the jump in 
the quantity, which is a second integral of the compacton 
differential equation, to be zero. For the RH equation we 
have that 



■p+l 
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F+'-pifTf 



2p-2 



(12) 



p+l" p+l 

must be zero when we go from xq — e to xq + e. 

D^sc[{fr(.x)fP-'ix)]^^^0. (13) 

This is always satisfied if there is no infinite jump in the 
derivative of the function. 

Following the same strategy starting with the CSS 
Eq. (2) one obtains instead the following equation for 
the solitary wave: 

^/^-^/'-o(/')V^ = 0. (14) 

We find ioT I — p + 2 the solution is again of the form 

f = Acos^/P[f3{x ^ ct)] , -tt/2 < f5y <tt/2 . (15) 

and zero elsewhere. In particular for {l,p) = (3, 1) with 
the choice a = 1/2, we find 



/3 = 1/(2^3) , A = 3c . 



(16) 



We find that both the RH and CSS equations allow 
for anticompacton counterparts to the compacton solu- 
tions, where anticompacton solutions correspond to the 
transformation / — ^ — /. We notice that the amplitude 
is proportional to the velocity c and that the conserved 
quantities M,Q in the RH equation, see Eq. (4), and 
M,E in the CSS equation, see Eq. (5), are proportional 
to odd powers of c. Thus in the scattering problem of 
compacton-anticompacton, these four conserved quanti- 
ties are zero. By convention we will choose the com- 
pacton to be traveling to the right with velocity c, and 
the anticompacton to be traveling to the left with veloc- 
ity — c and opposite amplitude. 

Whereas the properties of the RH compactons have 
been studied numerically extensively in the past [7, 22, 
33-43], the properties of Lagrangian-based compactons 
have only been discussed in the case of CSS compactons 
using Fade approximants [ ] and in the case of a fifth- 
order generalization of the KdV equation using pseu- 
dospectral methods [2S]. 

Following an earlier indication that a possible blowup 
may occur following a compacton-anticompacton colli- 
sion [2(S], in this paper we compare the case of ^(2,2) 
and L{3, 1) compacton-anticompacton collisions. Explic- 
itly, the K{2, 2) equation reads 

Ut + 2uU^+6Ua;U^x+'2uUxxx^0, (17) 

whereas the L{3, 1) equation for q = 1/2 is 

Ut + UUx + 4:Ux Uxx +2u Uxxx = . (18) 

The functional form of the K(2^ 2) and L(3, 1) compacton 
solutions is the same up to numerical coefficients. Our 
numerical study is based on the Fade methods discussed 
in Ref. ' These methods were used recently to study 
the stability and scattering properties of i4r(2, 2) [ ] and 
CSS [44] compactons. 

Although previous studies did not find substantial dif- 
ferences in the behavior of compacton-compacton scat- 
tering between the solutions of the RH and CSS equa- 
tions, what we find here is quite startling. In the CSS 
equation compacton-anticompacton scattering leads to 
annihilation, whereas in the RH equation compacton- 
anticompacton scattering leads to blowup. 

The paper is organized as follows: In Sec. II we review 
the basics of our numerical strategy using Fade approx- 
imants. This approach was discussed in detail in our 
previous publications [43, 44]. We present the results of 
our simulations for the scattering of a compacton and 
its anticompacton counterpart in Sec. III. We summarize 
our main findings in Sec. IV. 



II. NUMERICAL APPROACH 

The fourth-order Fade approximants used in this study 
follow from a systematic derivation [43] of the Fade 
approximants method [ ^] for calculating derivatives of 
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smooth functions on a uniform grid. Our derivation con- 
tained as special cases the Pade approximants first intro- 
duced by Rus and Villatoro [40-42]. 

In this approach, one considers a smooth function 
u{x), defined on the interval x € [0, L], and discretizes 
this function on a uniform grid, Xm — mh, with m = 
0, 1, • • • , M, and h = L/M . Then, Pade approximants of 
order k of the derivatives of u(x) are defined as rational 
approximations of the form 
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+ 0(Aa;'=), 
+ 0(Ax'=) , 
-f ©(Aa;'^) , 
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where E represents the shift operator defined as 

E ' Urn — Um+k ■ 



(19) 
(20) 
(21) 
(22) 

(23) 



Even- and odd-order derivatives require approximants 
that are symmetric and antisymmetric in E, respectively. 

Because the RH and CSS equations involve only 
the first three- lowest derivatives of the function u{x), 
and the fourth-order derivatives enter only because we 
need to introduce the artificial viscosity term needed 
to handle shocks, the simplest fourth-order Pade dis- 
cretization can be derived involving only the grid points 
{xm,Xm±i,Xm±2}- In this case, we start with the finite- 
difference approximation corresponding to the third- 
order derivative, Um^\ which is only second-order accu- 
rate, and we focus on obtaining fourth- or higher-order 
accurate Pade approximants of Um , Mm "* , and itm" 



the subset of grid points, {xm,Xm±i,Xr 



on 

±2}- A suite of 
different fourth-order accurate approximation schemes is 
desirable in order to make sure that our numerical re- 
sults are independent of the peculiarities of a particular 
approximation scheme. Just like in our previous stud- 
ies [4.S, 44], we consider next several sets of approxi- 
mants that mix fourth-order accurate approximations for 
two of the derivatives Um , Vm'^ , and Wm""* , with a sixth- 
order accurate Pade approximant for the third one, de- 
noted as follows: (6,4,4) - this approximation scheme 
is an extension of the scheme introduced by Sanz-Serna 
et al. [.V'>, .H4] using a fourth-order Petrov-Galerkin finite- 
element method, and corresponds to a sixth-order accu- 
rate third-order derivative, Um"''; (4,6,4) - this scheme 
features a sixth-order accurate approximation for Um"*; 
and (4,4,6) - this scheme was introduced first by Rus 
and Villatoro [40, 41] and corresponds to a sixth-order 
accurate approximation for Um^\ We also consider a 
(4,4,4) scheme that was shown to minimize the extent 
of the radiation train in our previous K{2^ 2) and CSS 
studies. 



We introduce a uniform spatial grid in the interval x G 
[0, L] by defining the grid points Xm = mAx, with m = 
0, 1, • • • ,M and the grid spacing /S.x — L/M. Then, 
the fourth-order accurate Pade approximants discussed 
above are applied to discretize the spatial part of the 
K{p,p) equation 

Ut-CoU^+T] U^a:xx ^ + {^^) ^^a: ^ ^ ' (2^) 

and the CSS equation 



Ut- CQUx+rj Uxxxx 



/_ 1 V >x 
'Pi-'-'-l)x + ^,i-'^')xxx-^^ 



(25) 



where subscripts t and x indicate partial derivatives with 
respect to t and x, respectively. Here, u{x,t) is time 
evolved in a frame of reference moving with velocity cq 
and in the presence of an artificial dissipation (hyper- 
viscosity) term based on the fourth order spatial deriva- 
tive, rj d'^u/dx'^. The latter is needed to damp out explic- 
itly the numerical high-frequency dispersive errors intro- 
duced by the lack of smoothness at the edge of the dis- 
crete representation of the compacton (see e.g. discus- 
sion in Ref. ['2><]). In practice, the artificial dissipation 
must be chosen as small as possible to avoid significant 
changes in the numerical solution to the compacton prob- 
lem. Nonetheless, the hyperviscosity is responsible for 
the appearance of tails and compacton amplitude loss. 

Following Rus and Villatoro [2'')], we numerically dis- 
cretize the time-dependent part of Eqs. (24) and (25) by 
implementing the implicit midpoint rule in time. Corre- 
spondingly, we obtain for the K{p,p) equation. 
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and for the CSS equation 
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(27) 



Here we introduced the notations, 



.,n+l _ 



Um (tn) and 

Um{tn + At), to indicate evaluations at two dif- 
ferent moments of time. We assume that u„i{t) obeys 
periodic boundary conditions, i.e. UM{t) = uo(i). 
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III. RESULTS AND DISCUSSION 

In the following, we specialize to the case of the K{2, 2) 
equation {p ~ 2), 

Ut + 2uUx + Uxx + 2 u Uxxx — , (28) 
which has the exact compacton solution 

«K22(a;, t) = y cos^ [J i{x, t)\ , \^{x, i)| < 27r , (29) 
and the case of the i(3, 1) equation (p = 1, / = 3), 

Ut + UUx + -iUxUxx + "^-U Uxxx = , (30) 

which has the exact compacton solution {a = \) 

Wpii3(a;,i) = 3c cos^ ^(a;,t) , |C(a;, t)| < n/Stt , 

(31) 

where c is the compacton velocity and we introduced the 
notation 



^(x,t) 



xo - (c - co)i , 



(32) 



with xq being the position of the compacton maximum 
at i = 0. 

Using Eqs. (1) and (2), one can study the possibility 
that the above compactons, uk22 and Upua, have anti- 
compacton counterparts, i.e. we check that these equa- 
tions allow for solutions corresponding to the transfor- 
mation and indeed wc find that both equations allow for 
anticompacton counterparts traveling with a negative ve- 
locity. 

From Eqs. (29) and (31), we see that the functional 
form of the K{2, 2) and L(3, 1) compactons is the same 
up to numerical coefficients. Both compactons are sta- 
ble, as discussed in our previous studies [43, 44]. Individ- 
ual compacton propagation in the compacton comoving 
frame of reference with zero hyperviscosity indicates that 
the compactons propagate with the emission of forward 
and backward propagating radiation. The amplitude of 
the radiation train is suppressed by reducing the grid 
spacing. Ax, indicating that the radiation is a numeri- 
cal artifact due to the lack of smoothness at the com- 
pacton edges. For a given choice of the time step, Ai, 
and grid spacing. Ax, the radiation trains emitted in the 
K{2, 2) and L(3, 1) compactons propagation show simi- 
lar trends and are self-similar. These results are robust 
with respect to the choice of the fourth-order accurate 
Pade-approximant numerical scheme [ ]. 

To reinforce the idea that until now the behavior of 
the if (2, 2) and i(3, 1) compactons was deemed vir- 
tually identical, despite the higher nonlinearity of the 
L(3, 1) equation relative to the K{2, 2) equation and the 
fact that the CSS equation is derived from a conserved 
Hamiltonian, we remind the reader that we have stud- 
ied also the pairwise interaction of compactons with dif- 
ferent amplitudes (velocities). Again, the scattering of 
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FIG. 1. (Color online) Collision of two RH compactons (top 
panel) and two CSS compactons (bottom panel) with c\ = 1 
and C2 = 2. Simulations were performed in the comoving 
frame of reference of the first compacton, i.e. co = ci, using 
the (6,4,4) scheme and a hyperviscosity, 77 = 2 x 10~^. Here, 
we depict the early development of the ripple created in the 
collision process. The RH and CSS collisions are qualitatively 
very similar. 



i(3, 1) compactons is very similar to the scattering of 
K{2, 2) compactons. In both cases, the compactons re- 
emerge with the same coherent shape. The collision pro- 
cess gives rise to a ripple originating at the collision site. 
The main difference here is that in the case of the scat- 
tering of K{2, 2) compactons we observe a shock being 
formed when the ripple switches from positive to negative 
values. No evidence of shock formation accompanying 
the i(3, 1) collision process was observed. In turn, the 
simulation of the L(3, 1) compacton scattering requires a 
much smaller artificial viscosity to obtain numerical sta- 
bility than in the case of K{2, 2) compacton propagation. 
We discussed in detail the pairwise compacton collisions 
in our previous studies [4.S, 44]. In Fig. 1 we compare the 
early development of the ripple generated in the collision 
process between two RH compactons (29) (top panel) 
and two CSS compactons (31) (bottom panel). In either 
case, the compactons have velocities ci = 1 and C2 = 2. 
The simulations were performed in the comoving frame 
of reference of the first compacton, i.e. cq = ci, using the 
(6,4,4) scheme and a hyperviscosity, rj — 2 x 10^^. As 
advertised, the positive- and negative-amplitude parts of 
the ripple decay slowly into low-amplitude compactons 
and anticompactons, respectively. Qualitatively the RH 
and CSS collisions are very similar. The process of ripple 
decomposition is much slower in the CSS than in the RH 
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FIG. 2. (Color online) Simulation of a K{2,2) compacton- 
anticompacton collision in the comoving frame of reference 
of the compacton, co — ci — — C2 — 1. We used the (6,4,4) 
scheme and a "variable" hyperviscosity, 10"'' < t] < 2x 10~^. 



FIG. 3. (Color online) Simulation of a L(3, 1) compacton- 
anticompacton collision in the comoving frame of reference of 
the compacton, i.e. co = ci — C2 = 1, using the (6,4,4) scheme 
and a hyperviscosity, 77 = 2 x 10^*. 



case, because the amplitude of the ripple is larger in the 
CSS coUision. 

In Fig. 2 we illustrate the collision between a com- 
pacton and an anticompacton (solutions) of the RH equa- 
tion (l = p = 2). Here, the compacton and anticom- 
pacton counterparts travel with velocities Ci = — C2 = 1 
and the simulation is performed in the comoving frame 
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FIG. 4. (Color online) Comparison of two simulations of the 
K{2,2) compacton-anticompacton collision, at t = 4, corre- 
sponding to hyperviscosities 77 = 10"^ and rj = 10"^. 
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FIG. 5. (Color online) Comparison of two simulations of the 
7^(2,2) compacton-anticompacton collision, at t = 7.5, corre- 
sponding to hyperviscosities r] = 10~^ and = 8 x 10~^. 
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FIG. 6. (Color online) Comparison of two simulations of the 
p = 1, / = 3 CSS compacton-anticompacton collision, at t = 3, 
corresponding to hyperviscosities r\ — 10^* and 77 = 10~^. 
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of reference of the compacton, i.e. cq = c\. We used the 
(6,4,4) scheme and a "variable" hyperviscosity, 10~^ < 
7y < 2 X 10~^. The collision process is symmetric and 
we observe the development of a series of shocks at the 
ends of the compacton and anticompacton located avi^ay 
from the collision site. The maximum amplitude of these 
shocks increases as new shocks are being formed. Hence, 
the simulations require larger and larger hyperviscosity 
values in order to overcome the instabilities formed at 
the shock fronts and to reach the later times in the scat- 
tering dynamics. The usable duration of the simulation 
depends on the size of the hyperviscosity parameter, 77. 
In order to track the compacton-anticompacton collision 
for t < 10, we performed simulations for a series of hyper- 
viscosity values in the range 10"" < ?7 < 2 X 10"^ The 
snapshots displayed in Fig. 2 correspond to the simula- 
tion with the lowest achievable value of hyperviscosity. 



FIG. 7. (Color online) Comparison of two simulations of 
the p = 1, I = 3 CSS compacton-anticompacton collision, 
at t = 7.5, corresponding to hyperviscosities 77 = 10"'' and 
rj = 8 X 10"^. The bottom panel illustrates the sensitivity of 
our numerical results with respect to the Fade discretization 
scheme. 



By time t — 8, the "smaller" shocks have disappeared 
and we are left with only one large shock front. From this 
point on, the amplitude of the compacton-anticompacton 
overlap blows up, indicative of a singularity being formed 
in the numerical solution. This behavior is similar to the 
blowup noted in the study of a fifth-order generalization 
of the KdV equation [2-^]. The latter numerical results 
were obtained using simulations based on pseudospectral 
methods. 
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In Fig. 3 we illustrate the collision between a com- 
pacton and an anticompacton (solutions) of the CSS 
equation {p = 1, I = 3). Again, we study the case of 
a compacton and its anticompacton counterpart travel- 
ing with velocities ci = — C2 = 1 and the simulation 
is performed in the comoving frame of reference of the 
compacton, i.e. Cq = Ci. Because, no shocks develop 
during this collision process, the simulation is performed 
in the comoving frame of reference of the compacton, 
i.e. Co = ci, using the (6,4,4) scheme and a hyper- 
viscosity, rj — 2 X 10^''. We observe that in the case 
of this CSS compacton-anticompacton collision the over- 
lap leads to an increase in amplitude with the maximum 
reached around t = 4, followed by a dramatic collapse. 

Finally, we comment on the accuracy of the numeri- 
cal results presented here. Because all simulations were 
converged with respect to the grid spacing. Ax = 0.01, 
and the time step, At — 0.002, our main concerns are 
focused on the need for the inclusion of a hyperviscos- 
ity, rj. As stated already, the artificial dissipation term 
is needed in order to damp out explicitly the numerical 
high-frequency dispersive errors introduced by the lack 
of smoothness at the edge of the discrete representation 
of the compacton: the compacton is continuous at the 
endpoints of its compact support, but the derivatives of 
the compacton at the endpoints are not. Unfortunately, 
the hyperviscosity gives rise to tails and compacton am- 
plitude loss, and the hyperviscosity value must be chosen 
as small as possible in order to avoid significant changes 
in the numerical solution to the compacton problem. 

In Fig. 4, we compare the numerical results obtained 
in the case of the K{2, 2) compacton-anticompacton colli- 
sion at t = 4, for hyperviscosities ?/ = 10~^ and rj = 10"'', 
respectively. Similarly, in Fig. 5 we compare the results 
obtained at t — 7.5, for hyperviscosities rj = 10~^ and 
?7 = 8 X 1G~^, respectively. 

In Fig. 6, we compare the numerical results obtained in 
the case of the CSS compacton-anticompacton collision 
at t = 3, for hyperviscosities 77 = 10"'' and 77 = 10^^ 
respectively. Similarly, in Fig. 7 we compare the results 
obtained at t = 50, for hyperviscosities 77 = 10"'* and 
?7 = 8 X 10^*^, respectively. Furthermore, for this last 
case, we also illustrate the sensitivity of our numerical 
results with respect to the Fade discretization methods 
discussed in Ref. [43], in order to demonstrate that our 
results are robust with respect to the choice of numerical 
discretization scheme. 



IV. CONCLUSIONS 

In summary, in this paper we used the Fade ap- 
proximants developed earlier for discussing the stabil- 



ity and scattering properties of the RH and CSS com- 
pactons [43, 44], to study the collision process of com- 
pacton with their anticompacton counterpart solutions 
of the K{2, 2) and L(3, 1) equations. 

Previous analysis of the stability and pairwise scat- 
tering of the K{2, 2) and L{3, 1) compactons suggested 
that the properties of these compactons were very simi- 
lar, with the only obvious difference being that the rip- 
ple generated at the collision site developed shocks in the 
K{2, 2) case, whereas no evidence of shock formation ac- 
companying the i(3, 1) collision process was observed. 
In turn, the simulation of the L{3, 1) compacton scatter- 
ing requires a much smaller artificial viscosity to obtain 
numerical stability than in the case of K{2, 2) compacton 
propagation. 

Our simulations of the collision of compactons with 
their anticompacton counterparts demonstrate that this 
scattering process leads to very different outcomes in the 
case of the K{2, 2) relative to the L{3, 1) case. In the 
case of the K{2, 2) collision, the compacton and anti- 
compacton develop a series of shocks at the ends located 
away from the collision site. These shocks appear with 
increased shock-front amplitudes and lead to the ulti- 
mate blowup of the if (2, 2) compacton-anticompacton 
overlap. The compacton and anticompacton counterpart 
do not re-emerge from the collision process involving solu- 
tions of the K{2, 2) and L{3, 1) equations. In contrast to 
the K{2, 2) case, in the L{3, 1) compacton-anticompacton 
collision process, the overlap collapses in amplitude, in- 
dicative of the different nature of the singularities being 
formed in the numerical solution in the two cases. 

We note that possible lines of further inquiry may 
concern a study of the compacton-anticompacton col- 
lisions for L(3, 1) equation corresponding to a = 3/2, 
in which case the K{2, 2) and L(3, 1) compactons have 
the same width. Also, it may be interesting to in- 
vestigate the compacton-anticompacton collisions in the 
L(p + 2;p) equations for p < I, because the compacton- 
anticompacton collisions for the K{p;p) show blowup for 
1 < p < 3. 
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